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Abstract. We report experimental measurements on the rheology of a dry granular material under a weak 
level of vibration generated by sound injection. First, we measure the drag force exerted on a wire moving 
in the bulk. We show that when the driving vibration energy is increased, the effective rheology changes 
drastically: going from a non-linear dynamical friction behavior - weakly increasing with the velocity- up to 
a linear force- velocity regime. We present a simple heuristic model to account for the vanishing of the stress 
dynamical threshold at a finite vibration intensity and the onset of a linear force- velocity behavior. Second, 
we measure the drag force on spherical intruders when the dragging velocity, the vibration energy, and the 
diameters are varied. We evidence a so-called "geometrical hardening" effect for smaller size intruders and 
a logarithmic hardening effect for the velocity dependence. We show that this last effect is only weakly 
dependent on the vibration intensity. 

PACS. 45.70.-n Granular systems. - 05.45.Xt Synchronization; coupled oscillators - 45.70.Mg Granular 
flow: mixing, segregation and stratification 

1 Introduction memory effect or in the fluid phase similar apparent rheol- 
ogy. Recently several propositions have emerged with the 

A wide class of materials a priori as different as colloidal a i m Q f describing in a general and unified way this phc- 

suspensions, pastes, emulsions, foams or granular mate- nomcnology (see for example [3E] and refs. inside). For 

rials, show very similar behavior in the vicinity of the granular matter under shear, the rheology in the dense 

jamming transition [J. In all these systems one may ob- p i iase i s essentially characterized by a friction coefficient, 
serve close to yield many common features such as aging, 
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i.e. a threshold model whose dependence with the shear- 
ing rate has been studied extensively However, in the 
presence of vibrations much less is known and it is often 
assumed that the mean kinetic energy should play a role 
similar to the thermodynamic temperature as for molecu- 
lar systems. This is essentially the case for granular gases 
where this analogy can be carried very far. However, this 
simple vision is currently challenged for dense granular 
matter. First, for dense granular packing under moderate 
shearing rates, the velocity fluctuations were shown exper- 
imentally and numerically to play a marginal role in the 
rheology [2]. Second, in the vicinity of the jamming tran- 
sition as for many other glassy out-of equilibrium phases, 
novel relations between fluctuation and dissipation were 
postulated to account for the specific modes of relaxation 
[5] . Many theoretical works suggested to replace the notion 
of thermal temperature by the more general concept of 
"effective temperature(s)" describing the structural evo- 
lution between blocked configurations (for a recent review 
see ref. [6] and refs. inside). For vibrated granular matter, 
several attempts were made to identify such relations in 
numerical simulations [7]. But experimentally, for many 
practical and possibly fundamental reasons, they essen- 
tially failed [819] due to the difficultly to probe fluctua- 
tions apart from energy input. The issue of approaching 
jamming can be viewed either from the solid elastic prop- 
erties (numerically see [12] and experimentally see |13ll4j ) 
or from the fluid regime [lOlllj : but so far, no consensus 
has been reached on the fundamental constitutive prop- 



Rheology of a sonofluidized granular packing 

erties and the rheological behavior of a granular packing 
under vibration. 

In this report, we access the effective rheology of a 
weakly vibrated granular packing under a constant con- 
fining pressure by measuring the force needed to drag a 
wire in the bulk at a constant velocity. We observe and 
characterize the different mobility regimes. To rationalize 
this situation we propose a simple model that reproduces 
some of the features observed experimentally In a second 
series of experiments, we measure the dragging force when 
the intruder is a bead attached to the wire and we esti- 
mate the dragging contribution induced by the presence of 
the spherical intruder when we vary the diameter and the 
level of vibration. Note that in the absence of vibration, 
similar intruder mobility experiments have been carried 
either in 3D [8116] or in a 2D granular packing [17] . 

2 Experimental setup 

The rheomeric device that we use here to explore the ef- 
fective material rheology is sketched on figQ] It consists of 
a glass container closed at its bottom by seven piezoelec- 
tric transducers that allow to input in a controlled way, 
a weak level of vibration in the granular packing. It has 
been shown by Caballero et al. [TS] that even a weak level 
of vibrations may create significant constitutive modifica- 
tions in the packing: this is the sono-fluidization effect. 

The container horizontal section is a rectangle of length 
L = 19.00(±0.05)cm and width W = 2.30(±0.05)cm. 
It is filled with a bidisperse mixture of glass beads of 
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Fig. 1. Sketch of the experimental setup. The driving is at a 
constant velocity V, The masses are such that Mi > M\. In 
inset is a typical resistance forces F(x) as a function of the 
horizontal displacement x 

1.0(±0.1)mm and 1.5(±0.2)mm mixed in equal masses 
(i.e. a mean bead diameter of d — 1.2mm), with the glass 
density being p = 2.6 10 3 Kg/m 3 . This value of bidis- 
persity was chosen to avoid crystallization. In a prelimi- 
nary version of the experiment, we had noticed that using 
mono-disperse grains could induce a crystallization which 
would show up as crystalline patterns visible on the outer 
surface as well as multiple plateaus in the compaction 
curve. To us, that was indicating a competing dynamics 
associated with the growth of crystallites, something we 
wanted to avoid a priori. Note that we did not perform real 
systematic measurements on this quite interesting feature. 
The size and number ratio we use here are just the results 
on an empirical test that seemed satisfactory in the con- 
text of the present report. The beads are in direct contact 
with the transducers. Each piezoelectric transducer is ex- 
cited by a square tension at a frequency / = 400-ffz that 
is set out of phase by tt compared to that of its neigh- 
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Fig. 2. Agitation energy measured with the accelerometer 
buried in the bulk as a function of input signal frequency. For 
these measurements the maximum possible electrical tension 
was applied to the piezoelectrics. 

bors. The external vibration driving is set externally via 
the electrical tension applied to the piezo transducers. We 
benchmarked the bulk agitation by monitoring the rms ac- 
celeration of a buried accelerometer that indicates a value 
7 = 0.35g (<? is the gravity acceleration) for the highest 
driving voltage. The agitation energy as a function of in- 
put signal frequency was measured with the accelerometer 
buried in the bulk. The results are shown in figure^ Note 
that the peak vibration intensity occurs around 400 Hz, 
which is the frequency at which the experiments were per- 
formed. On the packing top, we placed a thin metallic lid 
to probe via an induction detector, the surface position. 
The pressure applied by the lid corresponds to 5 granular 
layers: Pq = 5pgd = 150 Pa. 

This weak energy input in the granular material does 
not induce any motion that can be perceived by eye on 
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our current time scale, even if the accclcrometer shows 
that the system is being agitated. However, if we follow 
for a very long time scale the particle motion, we indeed 
observe substantial reorganization in the granular matrix 
[18]. Note importantly, that we have followed during 10 
days, at maximal input energy, a horizontal line of colored 
beads and we did not evidence any convection motion in 
the bulk. 

The wire is set in tension by two masses Mi < M-i. 
The masses values are chosen such that the mass differ- 
ence will always exceed the dragging force. Therefore, the 
mass Mi is always in contact with the force probe placed 
just below. The force probe which is a spring of stiffness 
k = 10 5 N/m, is fixed to a vertical translation stage driven 
downwards at a constant speed, V, by a brushless motor. 
We use two different types of intruders which are dragged 
horizontally and at a constant depth h below the surface. 
First, we use a metallic wire of diameter d = 0.1 mm, 
spanning the whole horizontal size of the cell, and second, 
we use the same wire but with a plastic bead of diameter 
D > d glued to it. The intruders are driven at constant 
velocities between V = 10~ 6 m/s and V = 5 10~ 3 m/s. 
Hence, the maximal value of the inertial numbers reached 

for these experiments is I max ~ Y ~ 10~ 3 . Hence, the 

y/P/p 

magnitude is much smaller than the values where shear in- 
duced hardening would be observed for sheared granular 
packing [2]. 

The experiments arc conducted using the following 
protocol. First, the packing and the intruder are vibrated 
at the maximal vibration energy during 4 hours. We checked 
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in few cases that using 10 hours would not lead to any de- 
tectable difference. At the end of this preparation phase 
a packing fraction of about <f> = 0.6 is reached. Then, 
the dragging experiment may start at a constant vibra- 
tion energy characterized by a kinetic energy fluctuation 
set externally by the voltage on the piezzo transducers. 
A typical force signal is presented in inset of figure {!]). 
Note that at the beginning of the pulling we observe an 
overshoot of the resisting force, but after few millimeters 
of displacement a steady-state is reached. 

3 Drag force on the wire 

In a first scries of experiments the effective rheology of 
the packing is probed by monitoring the drag force on a 
metallic wire. The wire is placed at a height h = 1.3cm 
below the surface. The average force necessary to drag 
the wire < F t >, is rescaled by the wire outer surface 
(St = 6 10 m) in order to obtain an average friction 
stress : a t = < Ft > /St- Then, the average friction stress 
is rescaled by the confining pressure P(h) = Pq + <f>p s gh 
to form an effective friction coefficient: /i e // = &t/P(h). 
This effective friction is displayed on fig. (T3]) as a function 
of the dragging velocity. We use two representations, a 
linear scale to display the lower velocity part of the mea- 
surements and in inset, a linear-log scale representing the 
whole range of measurements. The error bar corresponds 
to a r.m.s. deviation over 5 independent experiments. The 
data show several salient and robust features. First, at 
a constant driving velocity and for a small level of vi- 
bration, the pulling resistance decreases strongly with the 
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Fig. 3. Effective friction on the wire intruder /x e / / as a function 
of the pulling velocity V; no vibration (□), maximal vibration 
(•). In inset is the same graph but is a lin-log scale. The red 
line is the curve : y = x/8.5 10 -4 . 

vibration energy. Second, with or without vibration, the 
stress-strain relations arc of the hardening type: we in- 
deed observe an increase of the resistance to pulling with 
the dragging velocity. The third striking feature is that for 
the highest vibration energy and the lower driving veloc- 
ities, we reach a linear regime of the type: fi e ff ~ V/V* 
with V* ~ 8.510~ 4 m/s. Finally, for weak vibration en- 
ergy, the increase of friction with the driving velocity is 
weak, typically of the type [i = /io + /3hi(V/Vo), where Vq 
is a reference velocity. This behavior is also observed in 
the absence of vibration on about two decades in velocity. 
However, in this last case, we could not identify clearly 
the existence of a dynamical friction threshold as a well 
defined plateau. 

The inertial number being quite weak as noticed be- 
fore, the hardening effect observed here cannot be at- 
tributed to a shear induced hardening (at least in the 



usual sense as described in [5]). It is possible that the weak 
increase of friction with velocity comes from an intrinsic 
solid on solid friction property at the level of the gran- 
ular contacts. Note that such a behavior was evidenced 
for granular packing undergoing very slow shear against a 
wall [IS]. 

4 A drag force model 

Here we present a simple heuristic model suited to ra- 
tionalize the passage from a threshold rheology of the 
Coulomb type, to a linear force- velocity relation. The model 
shows explicitly that in the framework of solid friction 
coupled with elasticity, one can perfectly understand this 
change of rheology and the disappearance of a Coulomb 
threshold due to external vibration. Note that many so- 
phisticated variants of such an heuristic spring/friction 
model exist in the literature, essentially to account for 
the complex stick-slip dynamics of earth-quakes [19] . or 
to explain granular rheology fluctuations |17j . Under ex- 
ternal forcing, similar models were proposed in the engi- 
neering context to render the fatigue of structures |21I22| . 
In the presence of shearing vibration, the solid/solid fric- 
tion of a block on an inclined plane was investigated via a 
spring/friction model with memory [20j but the emphasis 
was rather on the stick-slip transition and the influence 
of mechanical vibrations on thermal aging. Note that it 
could be interesting in a future work, to bridge all these 
approaches and extend the present model to account for 
granular reorganization as well as contact thermal aging. 
A similar model was presented recently to account for the 
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Fig. 4. Rheology of the drag resistance model, coupling solid friction and elasticity under harmonic driving (see sketch on inset 
of On figure |3K is the rescaled average drag force < F > as a function of the rescaled pulling velocity V for an harmonic 
driving at Ud = 10 and for different vibration forcing amplitudes 7^ = 0.28(leftmost), 0.4, 0.7, 1.4, 2.8, and 7(rightmost). On 
figure |4jb) is the low velocity drag force threshold < Ftn > versus the vibration forcing amplitudes 7d for driving frequencies 
u>d — 1, 2, 10. The data collapse on the curve < F t h >= 1 — 7d (straight line). In inset, is displayed figure[4ja) in a lin-log scale 
to evidence the low velocity dynamical force threshold 

effect of a horizontally vibrating base on a spreading gran- force N results from the confinement pressure exerted by 
ular layer [10] . the medium on the intruder, and the friction coefficient 



A scheme of the model is sketched as an inset in figure 
0Ja). It consists of a linear spring of mass M in frictional 
contact with a block. On the block a constant normal 
force N is applied. The interaction between the spring and 
the block is through solid-solid friction characterized by a 
Coulomb friction coefficient /i, while the spring elasticity is 
characterized by a resonant angular frequency ujq . A force 
F(t) acts on the block such as to drive it at a constant 
velocity V. In analogy with the experiment, the block rep- 



accounts for the friction interactions. 

To render the influence of the piezo transducers ex- 
citing the internal elastic modes of the packing, we use 
a harmonic external driving force Fd(t) — fj,N Trfsin(wt) 
(see fig. HJ) j acting on the spring; the rescaled driving fre- 
quency is LUd = jj-, the non dimensional time is t = Lu^t 
and 7d is the vibration force amplitude rescaled by the 
friction force Fq = /iN [J 



resents the moving intruder while the spring accounts for 1 In the rest of the document, the tilde signs refer to rescaled 
the elastic response of the granular medium. The normal non dimensionalized quantities. 
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When the block and the spring stick together, we have characterized by an effective friction coefficient: 

x = V, where x = ^M^sl j s the non dimensional position of /t?\ , , t7 / ~ \ _1 

Veff = -jf = M {FJ = — I „ 2 _ j d I , (6) 
the spring and V = is the non dimensional velocity n X^d / 

of the block. In this case, the equation of motion is : which in terms of V s , eq. [5J can be written as 

F = S+F d , (1) W = ^I (7) 



as long as 
F = nN. 
When 



F 



< 1, with the forces F and Fd rescaled by 



= 1 there is a relative sliding velocity be- 
tween the spring and the block. In this case the spring 
equation of motion is: 



Equation (J7]) states that under strong external vibration 
and a very low velocity driving of the block, the effec- 
tive friction force that resists to motion increases linearly 
with the driving velocity. In this limit, the only impor- 
tant parameter that characterizes vibration is the r.m.s. 
x + F — Fd, (2) vibration velocity. That is quite consistent with what is 

observed experimentally under such conditions, as it is 
shown on figure [3J By numerically solving eq.[21 we display 
F sgn(V x). (3) on fig. |^a) the rescaled average drag force = as a 

Note that for simplicity, we do not include here any in- function of the rescaled velocity V for different vibration 
ternal dissipation mechanism that would damp the elastic amplitudes (labeled here in terms of the rescaled vibra- 



with Fd the external driving force, and 



resonance. 



For large vibrations, 



F 



« 



F A 



tion force amplitude 7^), at a constant driving frequency 
the block and the ^ = 10. We see that when the amplitude is increased, one 



spring always slide on each other (no sticking). In this obtains for ld ~ 1, a transition to a linear force velocity 

case, equation © reduces to the equation of a harmonic regime characterized by the limiting force- velocity relation 

oscillator driven by the external force %. It can then be 01 equation M The validity of limiting analytical solution 

easily solved. The velocity is then : (equation®) was verified explicitly by comparison with 



7 d cos uj d t, (4) 



the numerical solution. Another interesting feature is that 

iJJ<i x this model displays a transition between a threshold rhe- 

and the the r.m.s. spring velocity is : ology ftt bw vibrational driving) to a lincar force/velocity 

y~_ = ^j 2 ^} 1 ^ 2 = ^ U>d ~ d (5) rheology at higher vibration amplitudes. In the specific 

context of this model, the effective friction threshold to 

Considering that the block is moving at constant velocity 

trigger relative sliding can be written : 

V, one can derive from equations © and © the low ve- 
locity limiting value, V — ► 0, of the average friction force, ^thresh = A*(l — Id)- (8) 
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Thus the threshold vanishes for 7^ > 1 (sec figHJb)). Re- low as c = 40 m/s [T3], [2], where c « y/G/p, G being an 

lation ((SJ) can be simply understood by inspection of equa- effective shear stiffness. If we estimate the effective spring 

tion ([TJ) , where we see that the condition F < 1 is valid stiffness to be k « GLW/h, and the mass set into motion 

as long as x is smaller than x~m = 1 — 7^, which immc- as M = phWL, we get a resonance angular frequency 



diately yields the value for the effective friction threshold ujq = dkjM = c/h. The order of magnitude calculation 
(eq. [5]). This is in essence what seems to happen exper- of the resonance frequency fo ~ 500 Hz, is about what 
imcntally. It would be interesting in further experimen- we observe experimentally (see fig. [2]). 
tal work, to test more thoroughly the generic existence of 

this type of transition between a threshold fluid behavior 5 Drag force On an intruding bead 

and a linear fluid behavior for vibrated packing. At this 

In this second series of experiments, we measure the drag 

point, we do not seek to match exactly the experimental 

force exerted on an spherical intruder of diameter larger 

results, first because experimentally we do not have a clear 

than the grains. An essential difference with the previ- 

vision of which elastic modes are the most strongly cou- 

ous experiment is that now, the intruder motion induces 

pled to the wire friction. Also, because this model is not 

a back flow and thus, granular reorganizations in the bulk. 

suited to explain the slow (logarithmic) variations of the 

Pulling a spherical intruder at a constant velocity in the 

dragging force. This last behavior may be due to intrinsic 

bulk of a granular packing presents many practical difhcul- 

contact aging properties or unpinning dynamics in a dis- 

ties. Here, we try to override this problem by considering 

ordered environment. We just notice that the rms velocity 

the resisting force on an intruder with a bead attached to 

V s ~ 1.2 10 _3 m/s measured for the largest available vi- 

the metallic wire. The wire (with the bead attached) is put 

bration is of the order of the velocity V* = 8.5 10~ 4 m/s 

originally at an horizontal position ft, = 1.3 cm below the 

found for the effective friction relation in the linear regime 

surface and the resisting force on the system is measured. 

(fig. [3]) . This is quite consistent in magnitude with the pro- 

To extract the contribution of the bead Ff,, we also mea- 

posed theoretical picture. Finally, we just comment that 

sure the drag force on the wire only in an independent 

obtaining the elastic stiffness of the effective spring that is 

experiment at the same velocity and remove the corre- 

actually coupled to the wire, is a complicated issue. How- 

sponding mean resisting force from the drag force on the 

ever, we can provide an order of magnitude calculation. 

wire+bead intruder. Then, to obtain an effective friction 

Recent experimental measurements on surface waves in 

coefficient that characterizes the bead-only contribution, 

a granular packing under gravity have shown that under 

we divide Fb by the bead surface and by the confining 

such small confinement pressure, the propagation velocity 
for wavelengths of the order of one centimeter can be as 

subtraction procedure may be questionable as, in prin- 



pressure P : fj, e ff = F^/ttD 2 / ((f>pgh + P ). This simple 
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ciple, both contributions cannot be as easily decoupled. 
The presence of the bead may modify the overall friction 
stress on the wire with respect to the situation where the 
wire only is present. We nevertheless expect that this ef- 
fect will be small as the spatial extend of the perturbation 
due to the presence of the bead remains small compared 
to the wire length. The corresponding effective frictions 
are displayed on fig. El^a) for a intruding bead of diameter 
d = 2 mm. The error bars correspond to a typical r.m.s. 
dispersion over 3 independent experiments. 

We noticed that before reaching a steady state, typi- 
cally by pulling the intruder over a distance correspond- 
ing to a mean grain size d (this, almost independently of 
the intruder size), we go through a phase where an over- 
shoot of the dragging force is systematically evidenced. 
Also, similarly to the wire only case, the rheology is of the 
hardening type and the higher is the vibration, the lower 
is the resistance to pulling. 

In the absence of vibration we apparently reach, at 
low driving velocities, a threshold value smaller than its 
limiting value at higher velocities. Note that in previous 
dragging experiments performed in 2D, Geng et al. [17| ob- 
serve a very similar phenomenology. We estimate in this 
last case, the typical dragging velocities as being around 
10 _4 m/s. The confining pressure being of the order of the 
drag force in their experiment, they had an effective iner- 
tial number around 10~ 4 , which yields working conditions 
and a phenomenology consistent with the present report. 
But cither in their case or in ours, the reason for hardening 
effect observed at such low velocities is unclear. We can- 
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not dismiss, as in the case of the wire only, the existence 
of an intrinsic contribution coming from aging of granular 
contacts. However, note that in previous numerical works 
23 23] , it has been shown that material softness could also 
play a role as it leads to a hardening regime decoupled 
from the shear hardening regime for hard spheres. For the 
3D dragging experiment of a rod by Albert et al. [16) . no 
hardening effect was evidenced but the dragging velocities 
were essentially higher than our operating conditions. 

Interestingly and contrarily to the wire-only case, when 
the vibration is activated up to its maximal available value, 
we never could evidence a linear force- velocity regime, and 
this for any of the intruders we tested. We found essen- 
tially a slow logarithmic increase of the friction as dis- 
played on figGJi with D — 2mm on 4 decades in velocities. 

Now we vary the intruder size D. On fig. \Ejp and c 
we display, for different dragging velocities, the resisting 
forces without vibration (see fig©)) and for the maximal 
vibration (see fig. |Sfc) at different dragging velocities. The 
choice to represent directly the force instead of the effec- 
tive friction is motivated by the fact that the dragging 
force varies linearly with the intruding grains diameter D. 
Note that at the highest vibration level and for the highest 
drag velocities, we observe significant deviations from lin- 
earity. This linear Fb cx D relation is even more surprising 
if we choose to represent it as an effective friction coeffi- 
cient. Then, the drag force would be divided by the bead 
intruder surface (cx D 2 ) and we obtain a striking "geo- 
metrical hardening" effect for the smaller intruder sizes 
([leff oc -D -1 ). Interestingly, looking back at the data of 
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Fig. 5. Rheological properties of an intruding bead of diameter D. (a) Effective friction for an intruder of diameter diameter 
D = 2mm as a function of the driving velocity, for no vibration (□) and at the largest accessible vibration energy (•). (b) 
Dependence of the bead drag force contribution Fb as a function of size ratio D/d in the absence of vibration and (c) for the 
maximal vibration, (d) Estimation of the logarithmic hardening slope (3 es t as a function of the vibration kinetic energy for 
different diameters. 



Albert ct al. [16] , we indeed observe for the drag force on 
the cylinder, a crossover between a linear behavior (here 
corresponding to a constant effective friction) at large in- 
truder sizes, and a plateau corresponding to the "harden- 
ing" effect for smaller grain sizes. However, this effect is 
not present in the data of Geng et al. [T7] . 

Note that in the absence of vibration, the effective 
stress values are: fj, e ff(d = 2 mm) = 6 and /i e //(<i = 



6 mm) = 2. It is then straightforward to calculate the- 
oretically the friction force exerted on the surface of a 
sphere dragged in a medium for which the confining pres- 
sure is hydrostatic. If the local dynamical friction is //, 
the effective friction would then be: /x e // = ^^t. There- 
fore, the actual values for the effective friction coefficients 
found experimentally are quite large when compared to 
the usual values found for cohesionless granular packing. 
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This means that the effective friction should take into ac- 
count a flow processes in the bulk which involves many 
grains flowing around the intruder. It is clear that in this 
case, the lateral confinement could also play a role in the 
determination of the effective rheology, but still this does 
not explain why the effective friction would be enhanced 
for smaller grains sizes. 

Finally, we estimate the dependence of the logarith- 
mic hardening effect with the vibration energy. We seek 
to estimate the logarithmic slope of the effective rheology : 
[3 = dfi e ff /d\nV. In the context of solid on solid friction, 
this is a parameter quite sensitive to temperature which 
value increases significantly in the vicinity of the glassy 
transition [25] . Note that the whole protocol to determine 
a full rheological curve as in fig. [3j is quite cumbersome 
and lengthy. It took typically three months worth work. 
This explains why these results may look fragmented as we 
have not performed yet extensive and systematic measure- 
ments for a large range of driving velocities and vibration 
intensities. We rather choose to estimate the value of j3 
and the impact of vibration on its value. We actually per- 
formed two series of experiments at well separated driving 
velocities V = 3.3 10 _5 m/s and V = 5 10 _4 to/s. From 
those two points which average is determined over 5 in- 
dependent realizations, we estimate the logarithmic slope 
fi ea t- The results are displayed on figlSJi. We see that, in- 
deed, the increase of the vibration energy by a factor 5 
does not influence very much the velocity hardening ef- 
fect characterized by (3 es t- It is possible however from the 
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data, that in the absence of vibration the hardening effect 
would be smaller (typically a factor 2). 

6 Summary and conclusion 

In this report we presented experimental results demon- 
strating that even for a weak level of vibration, significant 
changes can be evidenced for the dynamical yield stress 
and the rheology of a granular packing. Two types of solid 
intruders pulled at constant velocity are used and we show 
that injecting sound waves in the bulk leads to important 
drag reductions. The resisting force generally increases 
with the pulling velocity, thus the rheology is of the hard- 
ening type. Interestingly, when a wire alone is driven very 
slowly, we observe a transition from a yield stress rheol- 
ogy characterized by an effective friction increasing slowly 
with velocity, to a linear force-velocity regime where the 
threshold vanishes for a finite energy input. The values 
obtained for the effective friction coefficient are consistent 
with a simple heuristic model coupling internal elasticity 
with solid friction. The logarithmic hardening regime ob- 
tained in the absence of vibration could be due either to 
vibrations generated internally by the bead friction mech- 
anism or, alternatively, by the intrinsic granular contact 
properties as evidenced by Ovarlez et al. [26j . Interestingly, 
when a wire alone is driven very slowly, we observe a tran- 
sition from a yield stress rheology characterized by an ef- 
fective friction increasing slowly with velocity, to a linear 
force-velocity regime where the threshold vanishes for a 
finite energy input. The values obtained for the effective 
friction coefficient are consistent with a simple heuristic 
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model coupling internal elasticity with solid friction. The 
logarithmic hardening regime obtained in the absence of 
vibration could be due cither to vibrations generated in- 
ternally by the bead friction mechanism or. alternatively, 
by intrinsic granular contact properties as evidenced by 
Ovarlez et al.[55]. We also measured the contributions to 
the drag force of spherical intruders of various sizes which 
motion induces a back-flow in the granular packing. In this 
case, logarithmic force- velocity relations were observed. At 
a constant drag velocity, we obtain a linear increase of the 
resisting force with the intruder size corresponding to a 
"geometrical hardening" effect that enhances the effective 
friction for smaller intruders. This is observed both in the 
absence of external vibration and at the maximal level 
of vibration. Finally, we estimate the logarithmic velocity 
hardening for different levels of vibration and found that 
it is only marginally influenced by the actual level of ki- 
netic energy. These two last effects remain a central issue 
to the understanding of granular packing rheology with 
and without vibration. 

We thank Jose Lanuza for help on the experimental device and 
Pr. Bruno Andreotti for discussion and critical reading of the 
manuscript. 
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